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Abstract 

i-j-". Deuteron Compton scattering is calculated to 0(Q 2 ) in pionless effective field theory using 

a dibaryon approach. The vector amplitude, which was not included in the previous pionless 
calculations, contributes to the cross section at 0(Q 2 ) and influences significantly the extracted 
values of nucleon electric polar izability at incident photon energy 49 MeV. We recommend future 
high precision deuteron compton scattering experiments being performed at 25-35 MeV photon 
energy where the nucleon polarizability effects are appreciable and the pionless effective field theory 
is most reliable. For example, a measurement at 30 MeV with a 3% error will constrain the isoscalar 
nucleon electric polarizability oo with a ~ 30% error. 
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I. INTRODUCTION 



Nucleon polarizabilities are fundamental properties of nucleons. They characterize how 
easily the nucleons deform under external electromagnetic fields. Forproton polarizabilities, 
tight constraints have been obtained by Compton scattering 0, 0, HI] and photoproduction 
sum rule (E For neutron polarizabilities, the extraction is more complicated and less 
satisfactory. Scattering neutrons off the Coulomb field of a heavy nucleus could in principle 
determine the neutron electric polarizability 0, IE 0]- However, the error is still quite 
large. By far the best methods known to measure the neutron polarizabilities involve photon 
scattering on deuteron targets. This introduces theoretical complications such as final-state 
interactions and two-body currents. 

Recently quasi-free Compton scattering (yd — > ynp) at 200 to 400 MeV photon energy 
was carried out at Mainz Q . Using the model developed in Ref. (§[ , tight constraints were re- 
ported on the neutron electric and magnetic polarizabilities: a n = 12.5±1.8(stat)ljjj;g(syst)± 
1.1 (model) and f3 n = 2.7 =f 1.8(stat)^(syst) =f 1.1 (model) in units of 10" 4 fm 3 (which will 
be used throughout this paper) jE]. The other approach, deuteron Compton scattering 
(yd — > yd), was carried out at Illinois [!(}, SAL (Saskatoon) [ll[ and Lund 12|| with photon 
energy from 49 to 95 MeV. Several theoretical calculations are performed using potential 
models [IE EE IE IE]. The latest calculation using the Lund data at 55 and 66 MeV gives the 
iso-scalar combinations of nucleon polarizabilities ao±/?o = 17.4 ±3. 7 and ao — /3q = 6.4 ±2.4 
which, when combined with the proton values, imply a n = 8.8 ± 2.4(total) ± 3.0(model) 
and (3 n = 6.5 q= 2.4(total) q= 3.0(model). 

An alternative method to analyze deuteron Compton scattering is nuclear effective field 
theory (EFT) (l7l EE] (see Ref. 0| for a recent review). The goal of the nuclear EFT 
program is to establish model independent, systematic and controlled expansions for few- 
body and eventually nuclear matter problems. For systems with characteristic momenta 
below the pion mass, theories wi th p ions integrated out are highly successful and applied to 
many observables (la.121. 22l 23l IE IE]. For systems with characteristic momenta above the 
pion mass, power counting is more complicated [IE]. But for practical purposes, the power 
counting developed by Weinberg [IE] could still be quite accurate despite its renormalization 
problem in the ^Sn channel. EFT calculations of deuteron Compton scattering were carried 
out in Refs. jEEE IE IE]- The latest extraction of nucleon polarizabilities using Weinberg's 
theory but with phenomenological wave functions yields «o = 8.9 ± 1.51 H a nd /3q = 2.2 ± 
i-5-a()i3]. In comparison, using a theory with pion integrated out, Ref . l2p e xtracts ao = 
8.4 ± 3.4 and (3 = 8.9 ± 4.3 from the 49 MeV data. Other calculations showed that 

data below 70 MeV give results consistent with the values predicted by leading-order chiral 
perturbation theory ao = 10 Po = 12 jib] ]. 

In the future, the High-Intensity Gamma Source (HIGS) at Duke University will be able 
to measure deuteron Compton scattering with high precision. Thus it is timely to ask 
what the best strategy is to improve the determination of nucleon polarizabilities. In this 
work, we focus on the low energy experiments — because the relevant theory is most well 
understood — and explore the sensitivity of ao and f3$ in future experiments. We follow the 
power counting in |28j, and present the complete unpolarized deuteron Compton scattering 
cross section to 0(Q 2 ) in the pionless dibaryon EFT j^, 24], 28||. A vector amplitude, which 
was not included in the previous pionless calculations was found to contribute at the same 
order as the nucleon polarizability contribution. The impact of the new contribution to the 
extraction of nucleon polarizabilities a and j3 is studied. 
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II. KINEMATICS 



The number of independent structures in Compton scattering amplitudes can be con- 
veniently analyzed in the helicity basis. In this basis, amplitudes are characterized by 
(h' x , h f 2 \hi, I12), where h^h^) is the helicity for particle i in the initial (final) state. Under 
parity transformation (P), a helicity amplitude transforms as 

P : (h[, h' 2 \h u h 2 ) -> (-h'i, -h' 2 \ - hi, -h 2 ) . (1) 

While under time reversal transformation (T), 

T : (h[, h' 2 \hi, h 2 ) -> (h u h 2 \h[, h' 2 ) . (2) 

It is clear that only the linear combinations of the helicity amplitudes that are invariant 
under P and T can contribute to Compton scattering. It is easy to see that there are 
2 ( J + 1) (2 J + 1) independent helicity amplitudes for a spin-J target. Thus there are two 
independent amplitudes for a spin-0 target, six for a spin-1/2 target, and twelve for a spin-1 
target, such as deuteron. The twelve amplitude structures for deuteron Compton scattering 
can be further classified as the scalar, vector and tensor amplitudes, S, V and T: 
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2M, 



N 



Se' d ■ e d + £ ijk Vie' d .e dk + T i5 (e' d .e dj + e di e' d . - -6^ ■ e d ^j 



(3) 



where e d (e' d ) is the initial(final) deuteron polarization, e is the proton charge and M N is the 
nucleon mass. By the same counting discussed above, S, V and T have two, four and six 
independent structures, respectively. We chose the following basis for these amplitudes |3l[ : 

S = f 1 e>*-e + f 2 s'*-s, (4) 
Vi = h(e'* x e)i + U(s'* x s)< + f 5 (s'* ■ ek, - e'* ■ skj) + / 6 (s'* ■ ek[ - e'* ■ sk<) , (5) 



Tij = / 7 < e. + / 8 sfs. + f 9 (V* • kk'^ + e ■ k!k t t'*) + f w (s'* • kkjs. + s • k'k. jS ;*) 

+/ u i (k,k, + k^k;.) e'* • e + / 12 i + KK) s'* • s - trace , (6) 



where e(e') is the polarization for the initial(final) photon and k(k') is the unit vector in the 
direction of the initial (final) photon momentum, s = k x e and s' = k X e'. 

For unpolarized Compton scattering, the amplitude squared is proportional to 

|A| 2 oc|5| 2 + ^|^| 2 + ^T^, (7) 

after summing over the initial and final deuteron polarizations. As we shall see in the 
following section, the \S\ 2 term starts to contribute at O(Q ) while a and (3q contribute at 
0(Q 2 ). The |Vi| 2 term also contributes at 0(Q 2 ) while the \Tij\ 2 term only contributes at 
0(Q 3 ). The \Vi\ 2 term should be included in a 0(Q 2 ) calculation to extract a Q and f3 
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III. UNPOLARIZED CROSS SECTION FROM EFT 



In the dibaryon formulation of the pion-less effective field theory , 24[ 28| , the nucleon 
field N and the 3 Si-channel di-baryon field tj are introduced. The effective lagrangian is 



£ = nUiD + 



D 5 



2M 
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iD + 
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AM 
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-y 



t)N T PjN + h.c. 



(8) 



where Pj = T 2 02O-i/v / 8 is the 3 Si two-nucleon projection operators and y is a coupling con- 
stant between the dibaryon and two-nucleon in the same channel. The covariant derivative 
is D = d + ieQA with Q — (1 + t 3 )/2 as the charge operator and A the photon vector 
potential. The nucleon-nucleon scattering amplitude is reproduced by the following choice 
of parameters 

8tt . 2/1 



V 



M 2 N r<? s ^ ' 



A 



M N r( is ^) Va (35l) 



(9) 



where a^ Sl ^ is the scattering length, r^ Sl ^ is the effective range, and /i is the renormaliza- 
tion scale in the power divergent subtraction scheme [3| which conserves gauge symmetry. 
Similarly, one can introduce the dibaryon field Si in the 1 5'o channel. 
The magnetic coupling in the lagrangian is 



2M N 
. e 



iVt (/i + fi lT3 ) o ■ BiV + 
L 2 



_ M N Vr( ls °M 3s ^ 



t]S 3 B j + h.c. 



M, 



N 



e^ttBjtk, 



(10) 



where /x — (A i p+A i n)/2 and \i\ = (fi p —fi n )/2 are the isoscalar and isovector nucleon magnetic 
moments in units of nuclear magneton, and B is the external magnetic field. Li has been 
determined by the rate of np — > dj [i^, |3(j. The measured cross section cr cxp = 334.2 ±0.5 
mb with incident neutron speed of 2200 m/s fixes L\ = —4.42 fm. L 2 is determined by the 
magnetic moment of the deuteron 



2 7 L 2 



1 _ 7r ( 3 50 



(11) 



in units of nuclear magneton, where 7 = \/MnB = 45.703 MeV with the deuteron binding 
energy B = 2.225 MeV. Fitting to the experimental value, one finds, L 2 = —0.03 fm. 

The relativistic correction of the above nucleon magnetic interaction is the "spin-orbit" 
interaction 



C 



so 
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CT-fDxE — ExD)N 
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where E is an external electric field. There are also polarizability contributions 
C pol = 27riVt ( ao + a t r 3 ) E 2 iV + 2nN^ % + fan) B 2 N 



2ixoia 



tUiE 2 + 



tU-B 2 



(13) 
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FIG. 1: Feynman diagrams for the scalar amplitude of deuteron Compton scattering. The thick 
initial and final state arrows represent the deuteron. The thick arrows in the middle of the diagrams 
are the dibaryon states in the 3 Si and 1 So channels. The open circles denote the electric current 
interactions, and the shaded circles the magnetic moment interactions. The solid squares are from 
electromagnetic counter-terms {L\ and L2). The open square is from the nucleon polarizabilities. 

where «o,i — ( a p=t« n )/2 is the isoscalar (isovector) nucleon electric polarizability, and 
Po,i — (Pp^Pn) /2 is the isoscalar (isovector) magnetic polarizability. a 4 and /? 4 are un- 
known two-body currents which will limit the precision of extractions of the above polariz- 
abilities in deuteron Compton scattering. 

In this version of EFT with pions integrated out, the ultraviolet cut-off scale A is typically 
set by the pion mass m n ~ 140 MeV, although the nucleon mass comes in from the nucleon 
propagator. Following ji^j, we will not distinguish the difference between Mjq and m n . 
There are some light scales denoted as Q in the system. The inverse scattering lengths 
1/aW and deuteron internal momentum 7 are of order Q while the effective ranges are 
counted as 1/Q. This counting allows the re-summation of the effective range contributions, 
which is a distinctive feature of the dibaryon approach. Note that in another approach, 
is counted as Q° [H, Hij], thus the effective range contributions are perturbative. 

The photon energy u) scales as Q 2 at energies comparable and below the deuteron binding 
(region I) and as Q at higher energies (region II) [22) . Since the nucleon polarizability contri- 
butions are proportional to u 2 , we will work in region II to accentuate their effects. Another 
scale in the problem is the relative momentum p between the nucleons in the intermediate 
states. Naively, p ~ y/M N u, thus one expects p ~ A when u> ~ 30 MeV. However, for the 
unpolarized deuteron Compton scattering, the pionless theory still converges well at u ~ 50 
MeV. This is because the main uncertainty from diagrams with nucleon-nucleon rescatter- 
ing in intermediate states are generally suppressed compared to the non-rescattering ones. 
Furthermore, the uncertainty is mainly from the P waves and higher, because the S-wave 
rescattering is described by the pionless EFT well in this energy range. 

Despite the successes of the pionless EFT, it is important that future experiments can be 
measured at lower energies such that the theory is unquestionably under control. We will 
explore the optimal energy range in the end of this section. 

Diagrams contributing to the scalar amplitude S (or equivalently fi and $2) to 0{Q 2 ) 
are shown in Figd The gauge-coupling diagrams in 1(a) and 1(b) are of order Q°, and of 
order Q 1 / 2 in 1(c). The magnetic diagrams in 1(d) is of order Q and of order Q 3 ^ 2 in 1(e), 
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while those in 1 (f )-l (i) are of order Q 2 , with the solid squares denoting magnetic two-body 
current L\ and L 2 . The open square in l(j) is the nucleon polarizability coupling which also 
contributes at Q 2 . 

The results of the diagrams are expressed in terms of the reduced amplitudes 



Si 



(l_ 7r ( 3 &)) 



Si 



(14) 



where i = 1 — 12. We write 



Si 

h 



gi(u) + gi{-u) + zh(uj) , 
g 2 (u) + g 2 (-u) - h(u) , 



(15) 



where 
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47 



2 7 ( 7 + 2V? 



TV 



3 ^7 + yj^f 2 — u Mjv — ie 



32 7T 7 a; Mjv «o -1 , V2 — 2 2 . 
H — ; _ — tan ( ; J 



e 2 y/2-2z 
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(16) 



and 

h(w) -- 
and where 

and 



l5u 2 M N A 



-4 7 ( 7 2 -cjMjv -ie) 2 + 87 6 + 207 4 c;M 7 v + 60 7 2 a; 2 M 
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-15 uj 3 M N 3 - 7 v/ 7 2 + Mjy (47 4 + 28 7 2 cjM iV + 39c; 2 M 7 v 1 



47 (7 - a/7 2 - w M N - ie^j ^ 



(17) 



327T7W M N f3 _ 1/V / 2 -2^, 
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jfc(u;) 
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uj L x M N + 2 I - 
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+ fc(w) 



N 



uj Mm — ie) H\ 



3M 



JV 



^-v/7 2 -^M w -z £ + ^(7 2 -^M w ) 



(18) 



(19) 



These results agree with those obtained in the previous calculations |28j, |29j. If the 
effective range is counted as order Q° instead of Q" 1 then the results of Ref. are 
reproduced. Also, the results of Ref. j2s| are reproduced up to recoil effects which can be 
safely neglected. 

The new contributions we calculate are shown in Fig. |21 Those are diagrams contributing 
to the vector amplitude V or, equivalently, the amplitudes /3,...,6- The solid circle denotes 
the spin-orbit coupling in Eq. (|12)1. The leading diagrams are 2(a)-2(c) which are of order 
Q. It is interesting to note that the relativistic correction (the spin-orbit interaction in 
diagram 2(a)) contributes at leading order. Diagram 2(d) is of order Q 3//2 , while diagrams 
2(e)-2(h) are of order Q 2 . Naively, if we want to calculate the cross section to order Q 2 , 
these orders Q 3 ^ 2 and Q 2 diagrams are not needed according to Eq. (JJJ). We still choose 
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FIG. 2: Feynman diagrams for the vector amplitude of deuteron Compton scattering. The solid 
circles denote the spin-orbit intercation in Eq. (|12|) . Other symbols are the same as described in 
Fig. 1. 

to include these contributions because there is an enhancement factor /i 2 ~ 5.5 ~ 1/Q in 
them. The diagrams give 

fa = 03 (w) - flf3(-w) + 2g(w) , 
/ 4 = flf 4 (o;) - g^-u) , 
7s = 0, 
/ 6 = -q(u) , 



where 



03 



27(2^0 + 2^-1) 



M. 



N 



tan 



(20) 
(21) 



V2-2^ " " v 4 7 
-2 7 5 + 2 7 3 y/j 2 + w Mjy - ie (27 - V7 2 - w Mjy - ie) 



(/jp + /ii) 

3 M w 3 vW^f 

+7wM w ^-77 2 + 37v / 7 2 + ^M w + 2^/7 4 -cj 2 M2 r -ie^ 



and 



+u 2 M N 2 (-5 7 + 3 vV + w Mjv 
3 

£ 4 (w) = - (/x 2 + /i 2 ) 



+ 



N 



Ml 



\J^ 2 — ijj M N — is 



(22) 
(23) 



with k(u) as given in Eq. (fTT)|) . 

As noted above, the tensor amplitude T only contributes to the cross section at (Q 3 ). At 
the same order, unknown two-body currents «4 and /? 4 defined in Eq. (fT3*|) also contribute. 
Thus we just calculate the unpolarized Compton scattering to order Q 2 . This will allow us 
to extract the nucleon electric polarizability ao with ~ 30% theoretical uncertainty. 

Combining the above results, the differential cross section for unpolarized Compton scat- 
tering is 



da 



a 



6(cu + 2M N ) 



; Re [3(1 + z 2 ){nh + f*f 2 ) + Uzflh 



+2(3 - z 2 )U;h + flU) + 8z/ 3 7 4 + 4(1 + 3z 2 )(f;f 5 + fife) 
+8z(3 + z 2 )r b U + 8(1 + z 2 )(f;U + flh) + 16z(/ 3 7 5 + fife)} 



(24) 
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9 

FIG. 3: Differential cross section for jd — * jd with incident photon energy u = 49 MeV (upper plot) 
and u =55 MeV (lower plot). The dashed curves are 0(Q 2 ) results with no nucleon polarizabilities, 
®o = 0o = 0. The solid curves are results using polarizabilities «o = Wfto = 1-2 X 1CT 3 fm 3 , 
calculated in leading order chiral perturbation theory j^. The dotted curves are the vector 
amplitudes contributions alone. The 49 MeV data is from Ref. ^J], while the 55 MeV data is from 
Ref. El- 



where a = 1/137 is the fine structure coupling constant and z = cos 8 = k • k . 

In Fig. 01 we show the differential cross section for u = 49 MeV (upper plot) and uo =55 
MeV (lower plot). The 49 MeV data is measured at Illinois while the 55 MeV data is actually 
a superposition of the 54.6, 54.9 and 55.9 MeV data from the Lund experiment. In both 
plots, the dashed curves are 0{Q 2 ) results with no nucleon polarizabilities, «o = Po = 0. 
The solid curves are results using nucleon polarizabilities calculated in leading order chiral 
perturbation theory (ChPT) |2o| . 

«o = 10/3o = ,J 9 }f 2 = 1-2 x lCr 3 fm 3 . (25) 

The dotted curves are the contributions of the vector amplitudes which were not included 
in the previous calculations 27, 2^]. These effects are ~20% and are as large as the 



nucleon polarizability contributions in the backward angles. 

In Fig. HJwe plot the cross section using the counting employed in Ref. [1^ and the value 
of Eq. (J25j) , We set 1/a^ ~ uj ~ Q, r« ~ Q° (such that the effective range contributions 
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FIG. 4: Differential cross section for — > "fd with incident photon energy to = 49 MeV (left plot) 
and u) =55 MeV (right plot). The leading-order ChPT values «o = 10/?o = 1-2 X 10~ 3 fm 3 are used 
and the effective range parameter is treated as order Q°. 
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FIG. 5: Differential cross section for jd — > 'yd with incident photon energy to = 49 MeV (left plot) 
and to =55 MeV (right plot). The best fit value «o = 12.3 and /3o = 5.0 (in units of fm~ 4 ) are used 
and all data points are assumed to be independent in the fitting. 



are treated perturbatively) and expand our amplitudes in Q. The leading order ChPT value 
of nucleon polarizabilities also gives a good description of data. Setting the vector amplitude 
to be zero, our result coincides with that of Fig. 2 of Ref . jl^ . 

Assuming the data measured at u =49, 54.6, 54.9 and 55.9 MeV are independent, a \ 2 fit 
gives ao = 12.3 ± 1.4 and (3$ = 5.0 ±1.6 (note that theoretical uncertainty is not included) 
in the dibaryon formalism. The best fit curves are plotted in Fig. |SJ In the counting that 
treats the effective range perturbatively, the x 2 fit gives a = 14.2 ±2.1 and j3 = 9.3 ± 2.5 
(theoretical uncertainty is not included) with very similar curves as those shown in Fig. 
El The difference between these two counting schemes is of higher order, thus it appears 
that while the fit values of a are quite close to the one predicted by leading order ChPT, 
the error on /3q is quite large using the current data. These results can be compared with 
«o = 8.4 ± 3.4 and /?o — 8.9 ± 4.3 from Ref. j2^|, where the contribution of the vector 
amplitude has been neglected. Thus the vector amplitude has significant effect on a reliable 
extraction of «o. 

Future high energy, high precision experiments can take advantage of the theoretically 
clean pionless EFT by performing the measurement at lower energies. However, the nucleon 
polarizability contributions become less important at lower energies. To study the best en- 
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FIG. 6: Nucleon polar izability effects on Compton scatterong cross section with ao = 10Po = 12 
(in units of 10 -4 fm 3 ). The curves, from top to bottom, are for photon energy u> = 20, 30, 40 and 
50 MeV, respectively. 

ergy range for future experiments, we show the nucleon polarizability effects in the Compton 
scattering cross section for «o = 10/?o — 12 (the leading order ChPT value) in Fig. El The 
curves, from top to bottom, are for photon energy uj = 20, 30, 40 and 50 MeV, respectively. 
Given the unknown 0(Q 5 ^ 2 ) effect is about 5% for a series with expansion parameter ~ 1/3 
and the theory converges better at <30 MeV, we recommend the best energy range to per- 
form the deuteron Compton measurement is 25-35 MeV. From this figure, we see that a 
measurement at 30 MeV with 3% error will constrain ao with 3 x 10™ 4 fm 3 experimental 
error, which is comparable to the expected theoretical error. 

IV. CONCLUSION 

We have computed the unpolarized Compton scattering in dibaryon effective field theory 
at 0(Q 2 ). The vector amplitude, which contributes to the cross section at 0(Q 2 ) but was not 
included in previous pionless EFT calculations, was found to affect the extraction of nucleon 
electric polarizability by more than 50% at 49 MeV. We recommend future high precision 
deuteron Compton scattering experiments be measured at 25-35 MeV photon energy for 
appreciable nucleon polarizability effects and controllable theoretical higher-order effects. 
More specifically, a measurement at 30 MeV with 3% error will constrain a with 3 x 10~ 4 
fm 3 experimental error, which is comparable to the expected theoretical error. 

This work was supported by the U. S. Department of Energy via grant DE-FG02-93ER- 
40762 and by the National Science Council of Taiwan, ROC. JWC thanks Paulo Bedaque for 
organizing the Summer of Lattice Workshop 2004 at Lawrence Berkeley Laboratory where 
part of this research was completed. 
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